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Abstract 

We have developed a field theory formalism to calculate in-medium properties of hadrons 
within a unified approach that exploits a single Lagrangian to describe simultaneously both mat- 
ter background and meson fluctuations. In this paper we discuss the consequences on physical 
observables of a possible phase transition of hadronic matter taking place in the chiral limit. 
We pay special attention to the pion velocity Ujr, which controls, through a dispersion relation, 
the pion propagation in the hadronic medium. The Vj^ is defined in terms of parameters related 
to the matrix element in matter of the axial- vector current, namely, the in- medium pion decay 
constants, ft and fs- Both of the pion decay constants change dramatically with density and 
even vanish in the chiral limit when chiral symmetry is restored, but the pion velocity does not 
go to zero, decreasing at most 10% over the whole density range studied. A possible pseudogap 
structure is indicated. 
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1 Introduction 



The initial motivation for building high energy heavy ion colliders was the quest for the quark 
gluon plasma. However, the data and the theoretical developments have shown that the phase 
diagram of the hadronic matter is far richer and more interesting than initially thought. At 
high temperature and/or density, hadrons are expected to possess properties that are very 
different from those at normal conditions. Chiral symmetry, which under normal conditions is 
spontaneously broken, is believed to be restored under extreme conditions. The value of the 
quark condensate {qq) of QCD is an order parameter of this symmetry, and is expected to drop 
as the temperature and/or density of hadronic matter is increased. Since the pion is the "litmus" 
indicator for spontaneously broken chiral symmetry, the various patterns in which the symmetry 
is realized in QCD, will be directly reflected in the properties of the pions in medium. 

The most fundamental quantities governing the dynamics of the pion are its mass and 
its decay decay constant fj^. These quantities have been the subject of previous studies in 
the formalism adopted here |2|. The property we are interested in is encoded in the pion 
dispersion relation in medium, which involves, besides the mass, the so-called pion velocity in 
medium, v-,^. This allows us to gain more insight into the real time properties of the system 
under extreme conditions and enable us to analyze how the phase transition from normal matter 
to deconfined QCD phase takes place bottom-up from the hadronic side. This issue has been 
addressed recently in heat bath IS [3] . 

At nonzero temperature and/or density, we have to take into account that the Lorentz 
symmetry is broken by the medium. In the dispersion relation for the pion modes (in the chiral 
limit) 

pI = vI\p\\ (1) 

the velocity v.,^ which is 1 in free-space must depart from 1. This may be studied reliably at 
least at low temperatures and at low densities via chiral perturbation theory [2]. The in-medium 
pion velocity can be expressed in terms of the time component of the pion decay constant, /* 
and the space component, 

(0|A0|7r''(p))i,.^edium = ^/*5«y, 

The conservation of the axial vector current leads to the dispersion relation (pQ) with the pion 
velocity given by 

vl = K/fl- (3) 

Other authors define two decay constants, ft and fs, different from those of eq. (|2j), 
through the effective Lagrangian, 

Ces = ^TridoU^doU) - ^Tr{diU^diU) + • • • , (4) 

where U is an SU(2)-valued chiral field whose phase describes the in-medium pion. In terms of 
these constants, the pion velocity is defined by 

Vtt = fs/ ft- (5) 

What we are interested in here is the pion velocity in dense medium as one approaches chiral 
restoration at the critical density pc- Before entering into this matter, we briefly review the 
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status of at high temperature which is more relevant to heavy-ion colhsions and the early 
Universe. There are two basically different predictions depending upon what degrees of freedom 
are taken into account, namely at T = T^, = in one version which might be referred to as 
"standard" and ~ 1 in another version which is non-standard. This stark difference is the 
motivation for our investigation of dense matter. 

Consider the scenario described by the two-flavor linear sigma model where the relevant 
degrees of freedom in heat bath are the pions vf and a scalar a. In this model, as the temperature 
approaches the critical temperature Tc, chiral symmetry becomes restored, the pions and the 
scalar join in an 0(4) multiplet, the phase transition belonging to the 0(4) universality class. 
This is the standard scenario largely accepted by the community in the field. It has been 
shown [3] that this scenario predicts that at T = Tc, the pion velocity must vanish. The 
argument is simple. Elevating the isovector axial chemical potential fiA to an U{1) gauge field 
- a powerful trick, it follows from the low-energy effective action of QCD that the isovector 
axial susceptibility (ASUS) XA is proportional to that figures in eq. Q where ft is the time 
component of the fully renormalized pion decay constant in heat bath. While one cannot reliably 
compute ft at low orders, one can use the information that the ASUS must equal the isovector 
vector susceptibility (VSUS) xv at the phase transition. Now one knows that xv\t=Tc 
from lattice measurements which implies that ft\T=Tc 7^ ™d that fg\'j'^j'^ = from a general 
consideration. Thus we have at T = Tc 

Vn = fs/ft = 0. (6) 

This is a simple and unambiguous prediction of the standard linear sigma model. 

The situation is markedly different when other light degrees of freedom are relevant. It has 
been shown ^5^ that when light vector mesons are in the picture as predicted by the hidden local 
symmetry theory with the "vector manifestation" a la Harada and Yamawaki [Zj (referred to as 
HLS/VM), the pion velocity is close to 1 at all temperatures near Tc- The crucial element that 
leads to this result is that the HLS theory with the vector mesons p and the pions coupled gauge 
invariantly and matched to QCD at the matching scale A^ has the vector manifestation fixed 
point at which the parametric mass of the vector meson Mp and gauge coupling g vanish. To 
this fixed point flows the system as one approaches the criticality (in temperature, in density or 
in the number of flavors). To one-loop order in a generalized chiral perturbation theory in the 
presence of light vectors ~ expected to be a good approximation near the critical point, the VM 
assures that XA = Xv at T = Tc, {f!^, f^) — > (0,0) (with /tt's defined in ((21)) as T — > Tc. As a 
consequence ^ 

~ 1. (7) 

We now turn to the main objective of this paper, which is to see whether there is a similar 
dramatic dependence of the pion velocity on the degrees of freedom also in dense matter. At 
present, although it has been established that hadronic matter at high density flows to the VM 
fixed point jH], describing dense matter in HLS/VM is found not to be easy because of coupled 
renormalization group equations in the renormalization scale and in density. Up to date, no 
conclusive result on physical variables (e.g., /tt, susceptibilities etc) has been obtained. We need 
therefore to resort to other formalisms. 

At low density, the standard chiral perturbation theory without vector degrees of freedom 
should be applicable. The in-medium pion decay constants have been computed to first order 

^The pion velocity in this scenario is not exactly equal to 1. There is a small deviation from 1 due to the 
fact that the heat bath violates Lorentz invariance by a small amount at the matching point Am- This has been 
computed [S]. 
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in the baryon number density p TTT. The results are = /7r(l — 0.26p/po) and = /^(l — 
1.23/9/po)) where po is normal nuclear matter density. Note that the space component is highly 
suppressed with respect to the time one. Therefore the in-medium pion velocity scales as 



The rapid drop of - and hence the pion velocity - suggests that leading order-order chiral 
perturbation theory cannot be trusted for density near that of normal nuclear matter. This 
observation calls for a careful re-examination of power counting. 

QCD sum rules have also been applied to calculate the in-medium pion decay constants jlll 
I12j . Unfortunately in the pseudoscalar axial- vector correlation function used, the dimension- 
5 condensate in the nucleon which cannot be calculated reliably is found to play the most 
important role in splitting the time and space components of the pion decay constant. This 
renders quantitative estimates difficult. For instance, a positive value of the condensate requires 
f^/fl > 1 and makes the pion mass tachyonic jlj, at odds with causality. In ref. [T2] . 
this problem was circumvented by taking the in-medium pion mass as an input parameter and 
obtained a negative condensate. It was found that the contribution of the intermediate A states 
is important in getting a sensible result. Using the input values for the in-medium mass in 
the range 139 MeV < m*^ < 159 MeV, the authors in [HI obtained /^//^ = 0.37 ~ 0.78 and 
= 0.63 ~ 0.79 and as a consequence the in-medium pion velocity in the range of 



In this paper, we study the in-medium pion velocity in the systematic field theory scheme 
developed in our recent works ^ |2j. There, the Skyrme picture is adopted to describe in a 
unified way both the pions and dense baryonic matter. A contact with HLS/VM theory of 
Harada and Yamawaki was suggested there. A static soliton solution having an FCC crystal 
structure [HI mUni CHI ^ 

is found to provide a classical description of dense baryonic matter. The pion is incorporated 
as fluctuating fields on top of this dense skyrmion matter. The interactions between the pions 
and matter modify the pion effective mass and decay constants. Since the background skyrmion 
matter breaks Lorentz invariance, the resulting space and time components of the pion decay 
constant become different. 

In the following section, we describe the Skyrme model approach used for investigating the 
properties of the pions in dense matter. A possible pseudogap structure is described in Section 3. 
In Sec. 4, we study the in-medium pion velocity by taking into account the background matter 
effect up to the second order in the interactions. We make a brief conclusion in Sec. 5. 

2 Model Lagrangian 

The aim of our series of investigations has been to develop the full dynamics of dense matter 
from a unique Lagrangian, which we have taken for simplicity of Skyrme type. Here, we use a 

slightly difFerent approach - but similar in spirit - to a crystal structure of dense matter was discussed by 
Diakonov and Mirlin fTTl. These authors studied the behavior of nuclear matter with density by implementing the 
nucleon-nucleon interaction by means of collective coordinate quantization of the skyrmion- skyrmion interaction 
and arrived at the conclusion that at twice or three time nuclear matter density a crystalline structure could arise. 
This calculation rends an additional support that a crystalline structure may indeed figure at high density, as has 
been shown within our approach 



vi = n/fi^i-p/po. 



(8) 



1/3 <v^<l. 



(9) 
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modified Skyrme model Lagrangian jl81 1191 QUI I21j . which incorporates in an effective way the 
scale anomaly of QCD in terms of a scalar dilaton field ^: 

^ = f (f ) Tr(a;.C/t9/^f/) + Id.xd'^x - ^ ((x//x)'(ln(x//x) " i) + i) (10) 

where U = exp{iT ■ tt/ fi^) £ SU (2) and x is the scalar field. The parameters in eq. (|lUj) are 
related to physical properties, f-,^ is the pion decay constant in free space, the corresponding 
quantity for the scalar field and its mass. For simplicity we will work with massless pions. 

The Lagrangian (jlUj) in the B = sector describes only mesons and yields the following pion 
and scalar field solutions in the vacuum, 

U = l, x = fx- (11) 

The fluctuations about this vacuum can be described by 

U^ = exp{iT-ip/f^), and x = fx + X- (12) 
Expanded up to second order in ip and %, the Lagrangian (|10|) takes the form 

£ = ^d^ifadf'y^a + '2^^^X^''X - h^lx^ + ■■■■ (13) 

This Lagrangian can also describe, a la Skyrme j^SJj baryons, i.e., the skyrmions, as topo- 
logical solitons if a proper stabilizing term, e.g., the Skyrme term 

= ^Tr [[U^d^U, U^d.Uf) (14) 

is added. The winding number associated with the homotopy of the static solution ?7o(r) is 
taken as the baryon number, 

B = i^e'^^ J SrTiiuld'UoUldWoUld^Uo). (15) 

The numerical results are sensitively dependent on the parameter values chosen for the 
Lagrangian. Given that the Lagrangian used is highly approximate, our results cannot be taken 
at their face values. However we believe that the qualitative structure of the results can be 
trusted. In this work, we will take the pion decay constant to be the empirical value, say, ~ 93 
MeV. As for the parameters associated with the dilaton, we have no direct empirical information: 
neither the mass nor its structure is experimentally identified. In the literature j24| I25j . they 
are adjusted so that the model fits finite nuclei as well as nuclear matter. The dilaton decay 
constant comes out to be ~ 240 MeV, while the its mass ranges widely from 0.5 GeV to 1.5 
GeV. We shall favor the lower masses, i.e. ~ 700 MeV, as proposed by several investigations 
|22| I26[ I27j although we shall report results for the wide range. 

The parameter e of the Skyrme term that figures crucially in stabilizing the skyrmion also 
plays an important role in setting the scale of dense matter. Unfortunately it is not obvious how 
to implement density dependence in this parameter. One may naively interpret the Skyrme term 
as arising when the p vector meson is considered to be much heavier than the scale involved and 
so is integrated out. This cannot be correct however because other mesons can contribute in the 

■^The role of the dilaton field in our model can be better understood as an interpolating field of the "soft gluon" 
degree of freedom that locks to chiral symmetry as described in e.g. |22| . The "hard or epoxy gluon" degree of 
freedom can be considered as integrated out with its efi'ects embedded in the coefficients of the Lagrangian. 
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heavy-mass limit and can destabilize the soliton. Furthermore, this term cannot be representing 
the effect of the p vector in dense medium since the p meson mass drops due to the vector 
manifestation |7J- On the other hand, the Skyrme term figures at extreme short distances as 
in the process of the proton decay a la Rubakov |2H]- Therefore it seems reasonable to think 
that the single Skyrme term subsumes a large number of massive degrees of freedom having the 
appropriate quantum numbers lying above the scale involved in the process. This implies that 
the constant e - representing short-distance physics - must be varying only slowly as a function of 
density. This is somewhat analogous to what happens in nuclear physics where the coefficients 
of zero-range counter terms are "universal" in that they are more or less independent of the 
density. See ^Hl where such a nuclear physics phenomenon is encountered. For our purpose 
here we simply take it as a constant. This assumption can be lifted in a model in which the 
heavy degrees of freedom are explicitly treated for soliton structure as in hidden gauge symmetry 
theory. In fact the explicit account of vector mesons in the skyrmion structure is found to be 
crucial for the bound state description of the pentaquark 0"*" in the soliton model [30] • Although 
we do not have a proof, we believe that the explicit presence of the vector mesons would not 
significantly modify the qualitative feature we are finding at and above nuclear matter density. 

In this work, we pick the value of e to fit, with the model, either the baryon mass spectrum j31j 
or the axial vector coupling constant qa |32j . The former fit leads to e ~ 5.45 (but in this case, 
/tt cannot have its empirical value, and comes out 64 MeV), while the latter fit produces e = 4.75 
(with at the empirical value). 

Two skyrmions are most attractive when they are relatively rotated by an angle vr about 
an axis perpendicular to the line joining their centers. Thus the lowest energy configuration of 
skyrmions for a given number density is in the crystal phase with all the neighboring skyrmions 
in the most attractive orientations. Moreover, the ground state of the scalar field x is shifted 
from the constant value to a static field configuration xo(^) in order to reduce the energy of 
the whole skyrmion-scalar system. 

3 Phase Transitions 
3.1 Chiral restoration 

We review and update the numerical results on skyrmion matter given in ref. pj. The aim 
here is to give the background for the calculation of the pion velocity discussed in the following 
section with the error committed in the scale factor in j2j corrected. 

In Fig.l we show the energy per skyrmion E/B of the skyrmion crystal as a function of the 
baryon number density p. The quantities in the chiral symmetry broken phase characterized by 
a non-vanishing (x) are drawn by filled symbols, while those in the chiral symmetry restored 
phase with {x) = are presented by unfilled symbols. In the numerical calculation, the size and 
energy of the skyrmion have been scaled in units of (c/tt)"^ and Gvr^/Tr/e, respectively. When the 
decay constants are fixed to /jr = 93 MeV and = 240 MeV, the numerical results depend only 
on the parameter = my^/ef-,^. We show the results for three different values of p^, namely 
1.25 (by triangles), 1.63 (by circles) and 2.26 (by squares), corresponding to rriy^ = 550, 720 and 
1000 MeV for = 93MeV and e = 4.75. The density is given in units of 0.015(e/^)^, where 
the constant 0.015 is chosen so that it equals the normal nuclear density 0.17(fm)~^ for /^r = 93 
MeV and e = 4.75. We should stress that the exact value of the density does not have a precise 
physical meaning, since the density scales with e, the Skyrme parameter, which as mentioned 
above, is not well defined. However we believe it to be a useful, albeit qualitative, guide to the 
magnitudes of the densities involved. 
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Figure 1: E/ B as a function of density. Numerical data for the {x) 7^ phase are presented by 
fihed symbols and those for the (x) = phase are presented by empty symbols. The density is 
in unit of 0.015(e/7r)^, which equals normal nuclear matter density po = 0.17 (fm)~^ when the 
conventional values, e = 4.75 and = 93 MeV, are used. E/ B is given in unit of Gvr^/^/e. 

As the density increases we reach a value at which the chiral symmetry restored phase has a 
lower energy than the chiral symmetry broken phase, i.e., a chiral phase transition takes place. 
Note that, even beyond the phase transition density, (x) 7^ is a local minimum. On the other 
hand, in the presence of the background matter, (x) = is always a local minimum whereas in 
free space, it was just an unstable extremum of V{x) ■ This phenomenon can be easily understood 
using the approximation introduced in ref. [2|, namely, taking x(^) ^ ^ constant xo- In that 
approximation, the energy per baryon E/B becomes 

E/B = {E2/B){xo/fxf + E,/B 

+^(xo) X [Volume occupied by a single skyrmion], 

where E2^a/ B are the contributions to the energy per baryon of the quadratic and quartic in the 
Lagrangian (|10|) and ^(x) is the potential energy density of the dilaton field. Now, it is evident 
that 

d'^{E/B) 



dixo/f: ^2 



= 2{E2/B), (17) 

xo=o 

and the system is quasi-stable at xo = as long as 2(E2/B) > 0. These local minima become 
true minima depending on the density, the former phase for low densities and the latter for 
densities above the phase transition. 

As an aside, we note that the numerical data for the (x) = phase can be fit by a simple 
curve 

E/B = ap^/^ + b/p, (18) 

where the first term comes from the Skyrme term contribution E4/B that scales as p^^^ and 
the second comes from the potential energy of the dilaton field. The coefficient a is almost 
independent of fi^. 
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The chiral phase transition takes place at the density, pc, where the phase characterized by 
(x) = s-iid that by (x) 7^ have the same E/B. The critical density depends sensitively on 
the parameter In Fig. 1, we show results for, /i^ = 1.25, 1.63 and 2.26 for which pc becomes 
1.32, 2.90 and 8.15, respectively. The figure drawn in a small box is the critical density pc as a 
function of /i^^. Note that it scales as The dashed line in the figure is not a fit but 0.67;U^ 
line just passing through the data point at p^ = 1.63. 

3.2 Pseudogap phase 

In Fig. 2 we show the average values of a = ^TrC/ and xl fx '^^^'^ space, which provide direct 
information on the structure of the crystal and can be used as order parameters of the phase 
transition. These data show that there is a "phase transition" before the chiral restoration 
transition, at the density at which the expectation value of a vanishes. Recalling that the 
expectation of a is proportional to the quark condensate, this may indicate a change in the 
chiral symmetry structure due to destructive interference of the phase of the quark condensate 
with a nonvanishing modulus. In the case of the massive pions, such a phase transition is not 
abrupt but continuous and this is the reason why we failed to notice it in our previous work [2]. 
We identify this phase transition with {a) = while (x) 7^ as of a pseudogap type. This 
phase persists in an intermediate density region, where the {x/ fx) does not vanish while {a) 
does jSH]- We denote the density at which the pseudogap phase transition appears as pp. Since 
the vanishing of a is possible because each skyrmion maps the chiral circle onto space in order 
to carry nontrivial baryon number, it may be that the pseudogap phase is an artifact of the 
skyrmion matter. We point out however that a similar pseudogap structure was proposed in hot 
QCD ISl. 

The two step phase transition is schematically illustrated in Fig. 2. 

(A) At low density {p < pp), matter slightly reduces the vacuum value of the dilaton field 
from that of the baryon free vacuum. This implies a shrinking of the radius of the chiral 
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Figure 3: /n^ dependence of two step phase transitions. 



circle by the same ratio. Since the skyrmion takes all the values on the chiral circle, the 
expectation value of a is not located on the circle but inside the circle. Skyrmion matter 
at this density is in the chiral symmetry broken phase. 

(B) At some intermediate densities (pj, < p < pc), the expectation value of a vanishes while 
that of the dilaton field is still nonzero. The skyrmion crystal is in a CC configuration 
made of half skyrmions localized at the points where a = ±1. Since the average value 
of the dilaton field does not vanish, the radius of the chiral circle is still finite. Here, 
(a) = does not mean that chiral symmetry is completely restored. We interpret this as 
a pseudogap phase. 

(C) At higher density (p > pc), the phase characterized by {x/ fx) — becomes energetically 
favorable. Then, the chiral circle, describing the fluctuating pion dynamics, shrinks to a 
point. 

In the case of massive pions, the chiral circle is tilted by the explicit symmetry breaking 
term. Thus, the exact half-skyrmion CC, that requires a symmetry between points with values 
a = +1 and a = —1 cannot be constructed and consequently the phase characterized by (a) = 
cannot exist at any density. However, {a) is always inside the chiral circle and its value drops 
much faster than that of {x/ fx)- Thus, if the pion mass is small enough, a pseudogap phase can 
appear in the model. 

In Fig. 3 we show the two step phase transition at various values of the parameter p,^. 
Contrary to pc, Pp depends only insensitively on p^. For small p^ ~ 1, pp ~ pc and its value 
is smaller than that for large p^^. As p^ increases, pc increases proportionally to p^ as was 
seen in Fig. 1, while pp saturates to 0.78. Since the density unit is (e/^)'^ and recalling that 
•\/2e/7r ~ mp, we can summarize the phenomena in terms of the parameters e and as, 

1. At density pp ~ 0.78 x 0.015(e/7r)^ ~ 0.004m^ , a pseudogap phase transition occurs. 

2. At density Pc ~ 0.67p^ x 0.015(e/7r)^ ~ O.Olm^, the chiral phase transition takes place. 
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with rriy ~ m 



4 In-Medium Pion Velocity 

The pion and scalar in dense medium can be described by fluctuations on skyrmion matter 
through the Ansatze 

U = VU^Uo{f)^, and x = Xo(r) + X- (19) 

Expanding up to the second order in the fluctuating fields, we obtain 

C = ^G''''{r)d^ipad''ip, + eabc^adiipbV^{r) 

+l9^x9^X - ^M(r)x' + P!,{f)xd,ipa. (20) 

This Lagrangian describes single particles moving under "local background potentials" provided 
by the static field configurations Uo{f) and Xo{^)- We have shortened the lengthy expressions 
for the local potentials introducing G"''(r), V^{r), M{r), Pa{r). Their explicit forms are shown 
in ref. j2]- Assuming background matter to be known, eq. (|2Uj) shows, when compared with 
eq. (fT^ . that a direct ip-x (not only <f'^-x) coupling might appear, which is absent in free space. 
Hereafter, we will drop the tildes for the fluctuating scalar field. 

In refs. ^H], by replacing the local potentials with their averages over space, we obtained 

Ces = \{G'''')d^Vad^^a + ^^^X^^X " \{M)x^ + ■■■. (21) 

We interpret this result as an effective Lagrangian for the in-medium pions and scalar with 
"intrinsic" parameters {G"'"') and (M), which correspond to the in-medium pion decay constant 
and the effective mass for the scalar respectively. In this naive approximation, Lorentz symmetry 
is still preserved. Thus, the pion velocity is 1 at this level. However, static local potentials will 
definitely break the Lorentz symmetry. In this work, we are going to evaluate the deviation of 
the pion velocity from 1 due to this background. 

Our strategy in performing the calculation is based on the following result. When local 
interactions with background matter are taken into account they lead to an effective Lagrangian 
for pion dynamics in the form of eq. © which defines the quantities necessary for the calculation. 
Let us be more specific. As far as the single pion is concerned, to lowest order in the fields, 
eq. Q with U = exp(iT • if/fn) becomes ^ 

J~-cfi = ^ (^y) ^oifadoifa - ^ (^yJ diipadiipa H . (22) 

Thus the procedure consists of calculating the propagators for the pion and the scalar fields from 
pn]) by solving for the background potentials. Then, matching these propagators with those 
obtained from the effective Lagrangian H22[). we have for the pion 

Zt Pq-Zs 

*If one substituted U = exp(ir ■ 0* / ft) with renormalized pion fields (^*, one would have 

1 , , 2 . . 
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where Z^^} = {ft,s/U?. 

Before proceeding to the calculation, we first obtain relations between the parameters, ft^s 
and /^'^, defined in Introduction. The effective Lagrangian @ yields the following axial vector 
current components 



(24) 



where the last equality is obtained by substituting U = exp(ir • (p/ f^^). When the expectation 
value of these axial vector current components is taken as in eq. @, we obtain the relations. 



h = fi, and /, = ^fifi. (25) 

Similar relations are derived in ref. for the real part of the complex pion decay constants. In 
what follows, the results will be given in terms of ft and fs- 

In order to obtain nontrivial results for the pion velocity in the medium we need to take 
into account the higher order effects of the background potentials. We shall proceed here in a 
perturbative scheme. To do so we decompose the Lagrangian into an unperturbed part, Cq, and 
an interaction part, Cj. The simplest way to do so is to take the free Lagrangian (fT3|l as Cq and 
the rest as Cj. Explicitly, 



_ 1 

2 



1 (M(f) - ml)x^ + Pi(f)x5i(/^a. (26) 



The presence of the local interaction potential G"'^{f) in the pion kinetic energy term makes 
the quantization process of the fluctuating pions nontrivial. The conjugate momenta of the pion 
fields ifa are given by 

n, = ^ = via. (27) 



We obtain 

where 

and 



n = no + ni, (28) 

Ho = i(9iV" + di^^di^n + hixx + diXdiX + mlx^), dSHK) 

ni = -Ci. (|2Hb) 



Note that the interaction Hamiltonian TCj is simply —Cj. The free propagators defined by Ti.Q 
and the interaction potentials appearing due to Tij are summarized in Fig. 4. In Fig. 4, C^^^l), 
for example, is the Fourier transformed of the local potential G"''(r): 

G^^{£) = T^ f d?re'^-''G''\r), (29) 

''box v'box 

where the integration is over a unit box of the crystal and Vbox is its volume. Due to the periodic 
structure of the crystal only discrete momentum values are possible. 

We show in Fig. 5 the diagrams used to evaluate the self-energy. Only the diagrams for S^^<^j^ 
appear. The symmetry structure of the skyrmion matter allows a nonvanishing self-energy only 
for a = b, One easily gets the corresponding diagrams for S^^^ or S^^. To first order, E^^^ is 
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(.PO,P) 



X X 

{Po,P) 
1 

Pq— — ITT'X 




+ieabcP- V%£) 




iq- P"(£) 




{po,p)\ y {qo,q) 



M{i) - ml 



Figure 4: Free propagators and interactions for the pion and the scalar fields in the presence of 
a background skyrmion matter. The energy-momentum conservation S functions are not shown. 




Figure 5: Diagrams used to evaluate the self-energy of the (pa propagation up to second order 
in the interaction. Here, b runs over 1, 2, 3 and the intermediate states run over all £ ^0. 
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nothing but TLi^l = 0). Since Ti.^^{0) = 0, no mixing between the fluctuating pions and the 
fluctuating scalar occurs. Thus, the pion propagator for ipa can be expressed as 

1 1 

(30) 



where we have used that the self energy to this order is S^^jy,^(po,p) = — p^(G""(0) — 1). The 
superscript "(1)" means that the quantities are evaluated to first order. Comparing these results 



with the propagator fl^ . we obtain ft = fs = fwy G^"-{0). Since G""(0) is nothing but the 
average of G'^"'{r) over the space, our calculation thus far reproduces the results of ref. |2j. To 
the same order, the self-energy of the scalar field is "^xx = — m^. Since it is constant, this 



self-energy modifies just the scalar mass from the free value to ^ M{0). 

Now, let us evaluate the second order diagrams shown in Fig. 5. Again, the symmetry of 
the background skyrmion matter does not allow nonvanishing off diagonal components such as 

(2) (2) 

S^j<^j,(a / b) and S^Va- We obtain 



, (31) 

Pq- {P + - 'm\ > 
= 0, (32) 

Here and in what follows, we economize on notation, \0^\^ = J2b 0^{0^)* (also for lower indices), 
and we sum over index c. Again, in spite of the '^ — x coupling term in the interaction Lagrangian 
H26|) , "^xva vanishes so that the pion propagator and the scalar propagator can be simply written 
as 

(34) 



respectively. 

The Pq and p dependence of S^^^ is not so simple as that of Ti^^\ Let's consider the small 
energy-momentum region. If we expand the self-energy Tl^p^^p^^ in powers of Pq and p, then the 
coefficients to second order in po and lead to 



f \ 2 1 /)2v(l+2) 



fJ 2 dp'^ 

f \ 2 1 ;52v(l+2) 

Js \ ]^ _|_ _ ^Va^a 



UJ 2 dpi 



= C^iO), (35) 
0,0 (. ^ * '"x J 



To obtain this result, we have used the observation that d'^Yi^^^'^^ /dpidpj is diagonal in the limit 
of po ^ and p^O. Note that the second order correction terms contribute only to fs/ fn and 
since they are all negative the pion velocity becomes v.,^ < 1. 
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In a similar way, the scalar mass gets a second order correction given by 



f = ml + 4i,-^)(0, 0) = M(0) - E f^- (37) 



The results of our calculation are shown in Figs. 6. Both of the pion decay constants change 
significantly as a function of density and vanish - in the chiral limit - when chiral symmetry 
is restored. However, the second-order contributions to the fs and /tt, which break Lorentz 
symmetry, turn out to be rather small, and thus their ratio, the pion velocity, stays v.,^ ^ 1. The 
lowest value found is ~ 0.9. Note, however, the drastic change in its behavior at two different 
phase transition densities. At the lower density, where the skyrmion matter is in the chiral 
symmetry broken phase, the pion velocity decreases and has the minimum at p = pp. If one 
worked only at low density in a perturbative scheme, one would conclude that the pion velocity 
decreases all the way to zero. However, the presence of the pseudogap phase transition changes 
this behavior. In the pseudogap phase, the pion velocity not only stops decreasing but also 
starts increasing with increasing density. In the chiral symmetry restored phase both ft and 
fs vanish. This result can be comparable to that of ref. p| found in the heat bath, where the 
pion velocity approaches 1 while both the spatial and temporal pion decay constants vanish at 
T = T,. 

As stressed in Sec. 2, because of the strong parameter dependence, e.g., on e and m^, 
the precise values at which these transitions occur may not be meaningful. What is relevant 
independently of the precise value of e and rriy. is the overall qualitative behavior of the pion 
velocity as a function of the density. A qualitative idea can be gained by the rough estimates: 
The pseudogap type phase transition occurs at 

Pp ~ 0.004m^(~ lApo with rup = 770 MeV), (38) 

with the minimum of the pion velocity at v-^ ~ 0.9 and the chiral phase transition takes place at 

Pc ~ 0.01m^(~ 2.9po with = 720MeV), (39) 

where the pion decay constants vanish. For the above estimates, we have used, for reference, 
the result of the analysis of the deeply bound states of pions in heavy nuclei j35j : 

{f:/f^)%=o.6po = 0.78 ± 0.05 (40) 
which is consistent with our naive estimates given above. 



5 Conclusion 

We have developed a formalism to calculate in-medium properties of hadrons within a unified 
approach which describes both matter and meson fluctuations with the same Lagrangian. In 
this paper, we have focused on the pion velocity, which controls through a dispersion relation, 
the pion propagation in the medium. For simplicity, we have chosen massless pions, i.e., the 
chiral limit and therefore the propagation is entirely controlled by the velocity. 

A novel aspect that was uncovered here in the skyrmion picture in the chiral limit is the 
presence of a psuedogap-like phase. This phase was not noticed in the previous work _2_ with 
massive pions. As density increases, the system first undergoes a phase transition to the pseudo- 
gap phase where the average value of a vanishes while the chiral circle preserves a finite radius. 
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Figure 6: In-medium pion decay constants and their ratio, the pion velocity. 



Roughly speaking, chiral symmetry is only partially restored in this phase. If we increase the 
density further, the radius of the chiral circles continuously decreases to zero and the chiral sym- 
metry is truly restored. The density range where the skyrmion matter is in such a pseudogap 
phase depends strongly on the scalar mass. 

Whether or not this pseudogap phase is an artifact of the Skyrme model at a classical level 
is not clear. The skyrmion comes out as a topological object in mapping the SU (2) space 
defined on the chiral circle into the space. As a natural consequence, the space occupied by the 
skyrmions would have a value for {a) which is smaller than the radius of the chiral circle. It 
would be very interesting to check on lattice whether such a pseudogap phase can be realized in 
QCD as suggested by Zarembo in heat bath 34 . 

The behavior of the pion decay constants, ft and fa-, drastically changes in the various phases 
of skyrmion matter. The constants vanish in the chiral symmetry restored phase in the chiral 
limit. The pion velocity obtained as the ratio of these two pion decay constants can deviate at 
most by 10% from its vacuum value 1, since the second order terms of the interaction between 
the pion and background matter which break Lorentz symmetry have a small effect in the decay 
constants. After achieving its minimum value at p = pp, this ratio increases asymptotically to 1 
up to the density where the pion velocity is well-defined. 

Our model calculation produces a complex phase transition scenario which seems to differ 
from the HLS/VM scenario in heat bath !5_, despite the fact that the pion velocity ultimately 
tends to 1 as in the HLS/VM theory modulo Lorentz symmetry breaking effect, and also from 
the linear sigma model in the 0(4) universality class 4 . 

Finally we should point out the short-comings in the calculation and further works to be 
done. First of all, ours is not a realistic calculation in that the model we have used may not 
be describing nuclear matter properly. Specifically we do not have the Fermi liquid structure 
of normal nuclear matter. Even so, we do nonetheless expect the qualitative features of our 
calculation to be more or less unaffected even after quantum fluctuations are introduced and the 
liquid structure is recovered. Next we must note that our calculation was made to second order 
in perturbation theory which may be justified at low densities but we do not have at present a 



14 



control of this approximation for larger densities. 

Other short-comings are that the kinetic energy contributions are missing from the calcula- 
tion PB] and that the subtlety in the interchange of the chiral limit and infinite volume limit to 
address phase transitions needs to be clarified. 

The immediate improvements to be made in the calculation are as follows. Apart from the 
quantum effects that would bring a liquid structure to the matter, it would be necessary to assure 
that the in-medium "vacuum" state is stable under mesonic fluctuations. For instance, we need 
to examine the eigen- modes and eigenvalues of eq.(20), focusing on zero-modes and negative 
modes if any. This is an interesting issue connected with the possible pion condensation etc. In 
this connection, we need to depart from the chiral limit. 

All these short-comings and improvements are under investigation and will be reported in a 
later publication. 
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